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ABSTRACT 

An analytical approach to the problem of replenishing task 
forces at sea is investigated utilizing queueing theory concepts 
and a random walk model. The type of operation considered 
consists of two different type supply ships servicing L combatants 
of different type. Replenishment times are assumed to be independ¬ 
ent random variables with service rates dependent upon both supply 
ship and combatant type. The total replenishment time of the task 
force is derived as the sum of specific individual replenishment 
times involving idle times of supply ships. With the additional 
assumption of exponentlality for the replenishment times, the 
probability density function of the total replenishment time is 
found in terms of the conditional Laplace transforms of specific 
type sequences of servicing and the probabilities that these 
sequences occur. Techniques for the evaluation of these quantities 
are discussed, however an explicit solution for the general case 
is not attained. As an example of these techniques, the.case of 
L = 2 is discussed. 
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I. Introduction 

The striking power, flexibility, and mobility of the naval 
task force make it a unique weapon of modern warfare. The ability 
of such task forces to realize their potential in carrying out 
their assigned missions is greatly dependent upon the resupplying 
at sea of such necessities as ammunition, fuel, and food. This 
resupplying of ships at sea is accomplished by underway replenish¬ 
ment and since it is of such importance, it is vital that the 
operation be professional and efficient. 

During an underway replenishment, the entire force is 
extremely vulnerable and therefore its safety is enhanced by 
reducing replenishment time to a minimum. Consistent with this 
goal of minimizing replenishment time, the Navy has long been 
involved in research and development of improved equipment and 
techniques for underway replenishments. Faster, more versatile 
supply ships, high speed transfer rigs, and vertical replenishment 
by helicopters are a few examples of the progress of this research. 
The fleet itself is constantly striving to improve its professional 
capabilities through operational training and exercises. Thus it 
appears that considerable effort has been extended in improving 
underway replenishments from an operational point of view. On the 
other hand, the replenishment operation also depends upon prior 
planning and efficient scheduling. This aspect of the operation, 
however, has not received adequate attention. If improvement in 
the efficiency of underway replenishments is to continue, 
increased emphasis must be placed on the analytical modeling and 
and study of the underway replenishment operation. 
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The basic idea of underway replenishment is a group of 
supply ships resupplying a group of combatant ships while underway 
at sea. The application of queueing theory appears obvious — 
n service channels (supply ships) providing services for (resupplying) 
m customers (combatants). The operation can also be viewed as a job 
scheduling problem — n machines (supply ships) performing operations 
on (resupplying) m jobs (combatants). Both concepts have been 
applied to the problem of analyzing the underway replenishment 
operation. 

Gordon and Copes [3] treat the operation as a job scheduling 
problem. General expressions are developed for the total waiting 
time for both combatants and supply ships and for the total time 
to complete the replenishment under the assumption of known, 
constant service times. Sequences to minimize the above quantities 
are found for the cases n = m = 2 and n = m = 3. However, since 
most of the factors describing the underway replenishment operation 
are not fixed, it seems intuitively unrealistic to treat service 
times as constants. Makino [7] examines a scheduling problem with 
service times treated as random variables. The generating function 
of the distribution of the total time is found for the cases of 
n = m - 2 and n = 3, and m = 2. This approach, as in the case of 
most scheduling problems in the literature (see for example [1], 

[10]), requires that all jobs must start on the first machine. 

This restriction severly limits this approach in modeling the 
underway replenishment operation. 

The literature yields an abundance of work on queueing theory. 
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Many of these, however, are based on the "Poisson-Exponential" 
operating assumptions which do not apply under-the restrictions of 
the underway replenishment process-(see for example[4], [6], [9]). 

The difficulty is that the operating assumptions provide a contin¬ 
uous stream of customers with Poisson arrival rate, while the 
replenishment operation prohibits further arrivals once the replen¬ 
ishment operation has started. McCullough [8] attempts to avoid 
the problem of the absence of further arrivals by allowing the 
customers already present to continuously repeat the service process. 
Steady state solutions for the mean number of units at a channel and 
the time required for a customer to complete all channels are then 
found by employing the technique of Koenigsberg [5 ] for cyclic 
queueing systems. The main objection-to this approach is that 
cycling does not actually occur and thus non-steady state results 
are desired. 

Waggoner" [1]] finds the Laplace transform of the non-steady 
state distribution of the total operation time utilizing queueing 
theory concepts and a random walk model. Service times are assumed 
to be independent exponentially distributed random, variables with 
service rates dependent upon the type server. This technique 
closely conforms to the restrictions of the underway replenishment 
operation and provides a useful analytical tool for the investiga¬ 
tion of this process. 

In actual replenishment operations, the type and number of 
both supply ships and combatants may vary. However, rarely will 
a task force be composed of only one type of ship. For the most 
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part, mission requirements and mutual-defense considerations 
dictate that task forces be composed of different type ships. A 
typical task force might consist of an aircraft carrier protected 
by two anti-submarine type destroyers and two guided missile type 
destroyers. It seems intuitive that the service rate for each type 
ship above would be different for a given type supply ship. There¬ 
fore the model developed.in [11] does not apply to the typical task 
force, since it assumes that service rates are the same for all 
combatants. The objective of this paper then is to obtain the 
non-steady state distribution of total replenishment time by 
extending the results of [11] under the more realistic condition 
that the service rate also depends upon the type of combatant 
involved. 
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II. The Underway Replenishment Operation 

An underway replenishment is an operation in which n supply, 
ships (n ^ 1) provide services for m combatants (m ^ 1) while underway 
at sea. Each of the supply ships can be a different type and can 
provide one or more services for a combatant. There can also be 
different types of combatants and usually each combatant requires the 
services of each type of supply ship. Service times are random 
variables with service rates dependent upon both the type of ship 
providing the services and the type of combatant being serviced. 

Prior to the operation, the order the combatants go alongside 
the supply ships is determined according to the requirements and 
capabilities of the particular ships involved. In essence, this 
forms queues behind each of the supply ships, although it is not 
necessary that these queues are physically filled with ships. On 
signal, the combatants rendezvous with the service force and those 
combatants assigned first to each of the supply ships proceed along-- 
side and simultaneously commence replenishing. Usually m > n and 
the remaining combatants are assigned to screening stations 'to 
provide defense for the force or to lifeguard stations, to provide 
rescue services for the ships actually engaged in replenishing. 
Although these ships are assigned supplementary tasks, they maintain 
their order in the imaginary lines behind each server. As the 
operation continues, each combatant moves through the different 
stations according to his particular assigned sequence. When a 
combatant has been serviced by all the supply ships, it leaves the 
queueing system — however, it does remain with the formation. When 
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all combatants have left the queueing system, the operation is 
completed and the combatants depart the area. 

As can be seen above, the following restrictions must be met 
by the model: 

1. A finite number of customers initially present and no new 
arrivals once the operation begins. 

2. Service times are random variables with service rate 
dependent upon both the server and the customer. 

3. Servers operating in parallel and in series, simultaneously; 
and 

4. A customer serviced by all servers departs the system. 
Finally, for simplicity, the problem is formulated for the case 

of two different type servers and a given number of different type 
combatants. 
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III. Formulation of the Problem 


The queueing system consists of two heterogeneous service 
stations operating simultaneously and independently. L different 
type customers are initially distributed in queues of size M and N 
at stations I and II and no additional customers arrive during the 
operation. 

Service times are assumed to be independent exponentially 

distributed random variables with service rates dependent upon both 

server and customer type. That is, the i^^ customer serviced by 

station I is serviced with service rate for i = 1,2,...,1 and the 
th 

j customer serviced by station II is serviced with service rate 
for j “ 1,2,...,L. 


Service begins simultaneously at both stations at time zero. 

Each station services all customers consecutively in predetermined 
order. A particular customer queues until serviced by one station 
and then immediately moves to the other queue, if any, for service. 
Transit times between stations are small compared to service times 
and are considered negligible. On completing service at this service 
station the customer leaves the system. For example, the second 
customer at station I queues until the first customer finishes 
service and is then serviced with service rate X^. Upon completion 
of this service, the customer immediately moves to station II and 
queues until his turn for service. The customer is then serviced 
with service rate departs the system. The initial set up 

and flow of customers is shown in Figure 1. 
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The operation is completed when all customers- have been serviced 
by both stations. Since the queueing system is finite, steady state 
does not exist. Therefore, the non-steady state distribution of the 
total time to complete the operation ia most significant. 


I II 



THE INITIAL SETUP 
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IV. The Model 


Waggoner [11 ] shows that the above queueing system can be 
represented by a two-dimensional random walk. The coordinates, 
[X(t),Y(t)] indicate the state of the system at any time t, where: 
X(t) = the number of customers served by station I in time t 

0 ^ X(t) 4 L 

Y(t) E the number of customers served by station II in time t 

0 £ Y(t) 4 L 

The random walk starts at the origin at time zero and moves by 
non-negative integer size steps until it reaches (L,L). Each 
sequence of steps represents a specific sequence of service opera¬ 
tions. A horizontal step to the right indicates that station I has 
finished servicing a customer and a vertical step upwards indicates 
a service completed by station II. 

A result of [111 is that not all paths from, the origin to (L,L) 
are allowed by the queueing system-.. This arises from the fact that 
the total number of customers serviced by a station in time t cannot 
be greater than the number of customers initially waiting for service 
from that station plus the number of customers which have completed 
service with the other station in time t. That is 


X(t) 4 M + Y(t) 
Y(t) ^ N + X(t) 


Therefore the random walk is bounded.by .the six lines whose 
equations are: 


X(t) 


0 

L 

M + Y(t) 
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Y(t) = \ L 

Ln + X(t) 

Definition 1 : 

The UPPER BOUNDARY consists of those points.along the line 
Y(t) = L with abscissa K, where M ^ K ^ L-1. 

Definition 2 : 

The RIGHT BOUNDARY consists of those points, along the line 
X(t) = L with ordinate K, where N ^ K ^ L-1. 

Definition 3 : 

The LOWER BOUNDARY consists of those points, along the line 
X.(t) == M + Y(t) , where 0 ^ Y(t) ^ N-1. 

Definition 4 : 

The LEFT BOUNDARY consists of those points along the line 
Y(t) = N + X(t) , where 0 ^ X(t) ^ M-1. 

Definition 5 : 

A PATH is a random walk which starts at the origin and moves 
by non-negative integer size steps to (L,L) and does not cross 
any boundaries. 

These boundaries and a typical path site shown in Figure 2. 
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FIGURE 2 


THE RANDOM WALK REPRESENTATION 


















V. Total Length of Time of a Path 

The total length of time of a path, T, is the sum of the times 
of the individual steps of the path. For each path there are L 
vertical and L horizontal.steps; however, not all of these steps 
contribute to T , since horizontal steps and vertical steps occur 
simultaneously. In order to determine those steps which contribute 
to T, some basic ideas are introduced. 

Definition 6 : 

th 

is the time it takes to accomplish the i horizontal step 
for i = 1,2,...,L. 

Definition 7 : 

th 

Yj is the time it takes to accomplish the j vertical step 
for j = 1,2,...,L. 

Lemma 1 : 

th 

The time of the i horizontal step, X^, begins at the termi- 
s t 

nation of the i-1 horizontal.step, unless the latter step ends 

s t 

on the lower boundary. If the i-1 horizontal step ends on the 

th 

lower boundary, then begins at the termination of the i-M 
vertical step. 

Proof : 

s t 

Consider any path from the origin to (L,L). If the i-1 

step does not terminate on the lower boundary, then the path is 

th 

permitted to continue horizontally and the i step begins at 
this time. A horizontal step of this-type is. shown in Figure 3. 


18 








If the i-1 horizontal step terminates on the lower boundary, 
then the path is not permitted to move horizontally and at least one 
vertical step is necessary before the path can again move horizon¬ 
tally. This “mandatory" vertical step begins on.the lower boundary 

at coordinates (i-l,i-M-l) and terminates at coordinates (i-l,i-M). 

th th 

Thus the i horizontal step begins at-the termination of the i-M 

vertical step. This is shown in Figure 4. 



LOWER 

BOUNDARY 


In any case, even though it is not necessary that the path 
continue to move horizontally, the fact that it can or cannot 
determines the starting point of the next horizontal step. 























Lemma 2 : 


The time of the vertical step, , begins at the termination 
s t 

of the j-1 vertical step, unless the latter step ends on the left 

s t 

boundary. If the j-1 vertical step ends - on the.left boundary, 
then Yj begins at the termination of the j-N^^ horizontal step. 

Proof: 

The proof is analogous to - the proof of Lemma 1. 

It is convenient to decompose the set of all possible paths into 
the following four mutually exclusive and exhaustive subsets: 

A: The set of all paths which pass through the point (L“1,L) 
and do.not pass through any points on.the lower boundary. 

B: The set of all paths which pass through the point (L-1,L) 

and pass through at .least one point on the lower boundary. 

C: The set of all paths which pass through the point (L,L-1) 
and do not pass through any points on the left boundary. 

D: The set of all paths which pass through the- point (L,L-1) 
and pass through at least one point on the left boundary. 

The above concepts are now applied.to the problem of determining 
the total length of time of any path. 

Theorem 1 : 

The total length of time of any path in set A is 


T, = I X. 

A i.i 1 
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Proof: 


Consider any path In set A* By definition such.a path must pass 
through the point (L-1,L) and therefore its last step must be hori¬ 
zontal and takes time X^. By Lemma 1 and the fact that the path 

does not pass through any points on the.lower boundary^ begins at 

s t s t 

the termination of the L-1 horizontal step* Similiarly the L-1 

horizontal step begins at the termination of the L-2^^ horizontal 

step and takes time X^^ Continuing, in this manner the process is 

traced back to the first horizontal step which begins at the origin 

and takes time X^^. Summing the times of these steps yields 

\ * • -+^2 + ^1 = ^ 3 ^ ^1 
which was to be shown. 

Theorem 2 : 

The total length of time of. any.,path in set B which last passes 
through the lower boundary at ordinate k is 
L k+1 

t/= Z X + I Y k = 0,1,2,...,N-1 

i=M4-k+l j=l ^ 

Proof: 

Consider any path in set B which last passes through the lower 

boundary at ordinate k. By definition auch a path must pass through 

the point (L-rl,L) and therefore the-last step must be horizontal and 

takes time Xj^. By Lemma 1, Xj^ begins ..at..the. termination of the 
s t 

L-1 horizontal step, unless the latter step ends on the lower 

s t 

boundary. Similiarly the L-1 horizontal step begins at the 
termination of the L-2^^ horizontal step and takes time ^ unless 
this latter step ends on the lower boundary♦ Continuing in this 
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s t 

manner the process-Is traced back to:the M+k+1 horizontal step. 

Since the M+k**^ horizontal-step terminates on the. lower boundary, 

s t 

by Lemma 1 the M4^k+1 horizontal..step^ begins at the termination 

s t 

of the k+1 vertical step and takes time 2 * Figure 5. 

Therefore the total time-the path-takes in-movlng-from (Mfk,k+1) 

> ' 

to (L,L) .is found by summing.the horizontal steps'pointed out above. 
This yields 

- W2 * Wi ■ 

i=M+k+l 



FIGURE 5 
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Since the path passes through at least one point on the lower boundary, 

it cannot pass through any points on the left boundary. Then by 

Lemma 2 the k+l^^ vertical step begins at the termination of the 

vertical step and takes time Similiarly the k^^ vertical step 

s t 

begins at the termination of the k-1 step and takes time Y^. 
Continuing in this manner the process is traced back to the first 
vertical step which begins at the origin and takes time (See 

Figure 5.) Thus the total time the path takes in moving from the 
origin to (>W-k,k+l) is found by summing the vertical steps pointed 
out above. This yields 


'k+1 




+ Y2 + Yi 


k+1 

I 

j=l 


The total time the path-takes in moving from the origin through the 
point (M+k,k) on the lower boundary to the point (L,L) then is found 
by combining the two results found above; that is 
L k+1 

T = E X + Z Y k - 0,1,2,...,M-1 

i=M+k+l j=l J 


which was to be shown. 

Theorem 3 : 

The total length of time for any path in set C is 


L 



Proof: 


The proof is analogous to-the proof of Theorem 1. 





Theorem 4 : 


The total length of time of any path in set D which last passes 
through the left boundary at abscissa h is 

h ^ 

= 2 : Y + E X h = 0.1,2.M-1 

^ l=N+bH ^3=1^ 


Proof : 

The proof is analogous to the proof of Theorem 2. 
Consider the following simple example: 

Example 1 : M = 5 N = 4 and therefore L = 9 

The total length of time of a typical path in set A is 


T 


A 


9 

Z X. 

T 1 


= X^ + 


X2-H 


X3 + 


X, + 




X. + 


X7 + Xg 


+ 


The total length of time of a typical path in set B with k = 2 is 


Z X. + 
i=8 ^ 


Z Y. 

j=i 




These paths are shown in Figure 6. 








FIGURE 6 
AN EXAMPLE 
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VI, Distribution of Total Replenishment Time 


The theorem of total probability is now employed to find the 
non-steady state distribution of the total replenishment time. As 
a preliminary step the sets B and D are decomposed as follows: 

: The subset of all paths in set B which last pass through 
the lower boundary at ordinate k for k = 0,1,2,,N-1; and 
: The subset of all paths in set D which last pass through 
the left boundary at abscissa h for h = 0,1,2,...,M-1. 
Theorem 5 : 

The Laplace transform of the probability density function of the 
total replenishment time is 


f,j,*(s) = E[e 


-sT 


N-1 _ k 

A|A] P(A) + Z ' E[e B |B ] P(B ) 
k=0 ^ ^ 


+ E[e"®^C|c] P(C) + 


-sT 

Z E[e D |Dj^] P(Dj^) 
h=0 


k h 

where T , T j T , and T are given above in Theorems 1, 2, 3, and 4 
ABC U 


respectively. 


Proof : 


By the use of Laplace transforms and the theorem of total 


probability it follows that 




= E[e''®''‘^A|A] P(A) + E[e ®'^b|b] P(B) 

+ E[ IC] P(C) + E[e”®^D|D] P(D) 
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N-1 M-1 

Since B = \J'B and D ~ U D and the fact that the B, 's 
k-0 ^ h=0 ^ 

and D, ’s are mutually exclusive, the Theorem follows, 
h 

The next step would'be to evaluate-the above conditional Laplace 
transforms and probabilities in order to find an analytical solution 
for the general problem. Unfortunately under the time constraints 
imposed in writing this,paper, the above result was-not explicitly 
attained. The structure of the problem and proposed direction of 
solution are given below. 

Event A means that the path ends horizontally and does not pass 
through any point on the lower boundary. In order,.for the path to 
avoid the first point on the lower boundary (M,0),,it is necessary 
that 

+ . . . + X^- > 

Similiarly in order to avoid the second point on the lower boundary, 
it is necessary that 

\ + • • • ■- Vl’ 

and finally to avoid.the last point on the lower boundary, it is 
necessary that 

Xi + . . . + X^_l >Yi-H . . . + 

Further, in order for the path to end horizontally, it is necessary 
that 

Xl+. . . + X^> Yi+. . .+Yj^ + Io+. . 

where is the idle time, that is the time the path sojourns 
at the point on the left boundary (k,N+k). That is 
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and 


N 

+ 


' [X. - E 

^il- 


' k+1 

N+k 

k-1 + 

MS x^ - 

E Y 

- E I ] for k = 1,2,...,M-1 

j=o J 

h“l 

= max(X,0) 

1=1 


where [X] 

Summarizing these conditions for event A yields 


f 


(1) A 


M+r 

E > 

r+1 

E 


for : 

i=l ^ 

j"! 


L 

L 


M-1 

E X. > 

E . 


+ E 

• -1 1 
i“l 


k=0 


for r = 0,1,...,N-1; and 


Event means that the path ends horizontally and last passes 
through the lower boundary at (Mfk,k). In order for the path to last 
pass through (lH-k,k) it is necessary that 


^1 + • • • + ^k + ''o 


"^k-l ^ "^1 ^k+1 


,th 


where is the idle time of the path at the i point on the lower 
boundary and is expressed similiarly to above, namely 

M + 

-■o - '’'l- 



i+1 

M+i 

1-1 


and = 

[ E Y - 

h=l ^ 

E X. - 
k=l ^ 

E 

j-0 


Now in order 

for the path 

to avoid 

the k+1®^ 


boundary (lfl-k+l,k+l), it is necessary that 


^k+1 


> Y 


k+2 
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Similiarly in order to avoid the k+2 point on the lower boundary, 
it is necessary that 

^k+l ^k+2 ^ \+2 \+3 


and finally to avoid the last point on the lower boundary it is 
necessary that 

^k+l ^ \+2 \ 

Further in order for the path to end horizontally, it is necessary 
that 

+ . . . + + . . . + + . . . + 

th 

however, since it is known that the path sojourns at the k point 
on the lower boundary we have 


k+1 

M+k 

k-1 


1 X. - 

i=l ^ 

I 

h=l ^ 

j=0 


Substituting this expression for Jj^ into the final condition above 
yields the condition 

^k+l ^ ^ \+2 


Summarizing the above conditions yields for k = 0,1,2,...,N-1 


( 2 ) 



M+k 

k-1 


k+1 

Z X 

+ z . 


Z 

i=l ^ 

h=0 

h 

j=l 

M+r 


r+1 


Z 


Z 

h 

i=M+k+l 

i 

j«k+2 

L 


L 


Z 

X > 

Z 

''j 

i=M+k+l 

1 

j=k+2 


Y^ , and 
for r = k+l,k+2. 


,N-1, and 
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By an analogous argument the events C and are described as follows 


(3) C = 


N+r 


r+l 


E 

i=l 

Y. > 

1 

Z 

j=l 

X. 

J 

L 


L 


E 

i=l 

Y. > 

1 

Z 

j=l 

X. 

3 


N-1 


k=0 


and for h = 0,1,2,...,M-1 


( 4 ) . <; 


N+h 

h-1 h+1 


Z Y 

i=l 

+ 2 I, < Z X, 

k-O '' J 

and 

N+r 

r+l 


Z 

i=N+h+l 

Y. > Z X. for r 

^ j=h+2 J 

= h+1, h+2 ,. 

L 

L 


Z 

i=N+h+l 

Y, > Z X, 

^ j=h+2 J 



,M-1; and 


In order to evaluate the probabilities of the above events, the 
analysis could proceed as follows 
N 

A . n ^ 

r-O 

where is the r^^ condition of (1) above. By complementation 

_ N _ 

A = yj A 

r=0 

“ ^0 ^ ^^1 ^ ^0^ ^ ^• • • ^^-1^ 

Since each of these terms is mutually exclusive 

P(A) - P(Aq) + P(A^ Aq) + . . . + P(^.Aq, . . . A^_^) 

and since P(A) = 1 - P(A) 

then P(A) = 1 - P(Aq) + P (A A^) + . . . + P(A^ A^, . . 
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Thus the problem of evaluating P(A) consists of evaluating the 
probabilities on the right hand side of the above equation. A similar 
procedure would be followed for the events C, and 

The next question is the evaluation of the conditional Laplace 
transforms for each of the events. For example, an evaluation of the 
conditional Laplace transform of T. given the set A could proceed as 
follows. By Theorem 1 


L 



for any path in set A. Therefore 

L 

E[e"®^A|A] = E[e'® \\k] 

= E[e'®^l|A] . . . E[e"®^L|A] 

and the problem reduces to the evaluation of these factors. A similar 
procedure would apply for the remaining conditional Laplace transforms. 
As an example of the above techniques the case for M = N = 1 is now 
discussed. 

Example 2 : M = N = 1 L=M+N = 2 

Solution: 

By Theorems 1, 2, 3, and 4 


Ta = Xi + 




T = Y + Y 
C ^2 
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The event A is described as follows 


Aq - (X^ > Yj) 

= (X^ + X^ > + Yj + I„) 

since it is known that all paths in A pass through the left boundary 
at (1,0), then and therefore 

A^ = (X 2 > Y^) 

Although the independence of events (A^) and (A^) could be used 
to compute the probability of A directly, the procedures outlined 
above will be used for illustrative purposes. 

A = (X^ > Y^, X^ > Y 2 ) 

and P(A) = 1 - P(Aq) - P(A^ A^) 

= 1 - P(X^ < Y^) - P(X2 < Y 2 ,X^ > Y^) 

= 1 - P(X^ < Y^) - P(X2 < Y 2 ) P(X^ > Y^) 
and utilizing (Al) of Appendix A to evaluate these probabilities yields 


P(A) = 1 - 


^1 


Al + 




or P (A) == 


(a 1 1^1^ ^^2 ^2^ 

By a similar procedure the probabilities of the remaining events are 
computed and are 


B = (Xi < Yi,X2 > Y^) , P(B) = 


Al ^2 


(Ai 


y 2" 
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C = (X^ < < Y 2 ) , P(C) = 


Xi A 2 


(Xi + y^) (X 2 + y2^ 


D = (X^ > Y^,X2 < Y 2 ) , P(D) = 


Ml I2 


(Xi + y^) (X 2 + y2) 


The conditional Laplace transforms are computed next. For event A 

E[e"®'^A|A] = (X^ > Y^,X 2 > Y 2 )] 

since and are independent this yields 

= E[e~®^l|X^ > Y^] Ete"®^2|X2 > Y 2 ] 


employing (A4) of Appendix A this yields 
E[e-®^A|A] = ^ ^ ^ 


^2 ^^2 *^2^ 


(Xj^ + s)(Xj^ + + s) (X2 + s)(X2 + y2 + s) 


The conditional Laplace transforms of the other events are obtained 
by a similar procedure and are 

Pt (Xt + y,) X„ (X- + y-) 


E[e‘®'^B|B] 


E[e ®^C|C] = 


E[e"®^D|D] 


(y^ + s) (X^ + y^ + s) 

(X 2 + s) (X 2 + y 2 + s) 

y^ (X^ + y^) 

y 2 2 ^ y 2 ^ 

(y^^ + s) (Xj^ + y^ + s) 

(jli2 + s) (X2 + y 2 + s) 

Xi (X^ + y^) 

y 2 (X2 + y 2 ) 


(Xj^ + s)(X^ + + s) (y 2 + s)(X 2 + y 2 + s) 













Combining these results with the probabilities for the events computed 
above yields 


f/(s) = 


Xi 


Pi U2 


H ^2 ^1 ^2 


(Ai+Ui+s) (A2+)i2'''s) (A^+s) (A^+s) 


(A^+y^+s)(A2+y2+s)(y^+s)(A2+S) 


^1 ^2 ^1 ^2 

(Ai+yf+s) (A2+y2+s) (y^^+s) (y2+s) 


^1 ^2 ^1 *^2 

(Ai+yf+s) (A2+y2+s) (A^^+s) (y2+s) 


^1 ^2 '^1 ^2 ( A 2 + y 2 + 2 s ) 

(A^+y^+s)(A2+y2+s)(A^+s)(A2+S)(y^+s)(y2+s) 


34 









BIBLIOGRAPHY 


1. Dudek, Richard A., and Teuton, Otis F. ^Development of M-Stage 

Decision Rule for Scheduling N Jobs through M Machines," 
Operations Research ^ 12, 1964: 471-479. 

2. Feller, William. An Introduction to Probability Theory and Its 

Applications ^ Vol. I. New York: John Wiley and Sons, Inc., 1957. 

3. Gordon, Bradley W., and Copes, Raymond F., III "An Investigation 

of Optimal Scheduling of Underway Replenishments." United States 
Naval Postgraduate School, Monterey, June, 1967. 

4. Gumbel, Harold. "Waiting Lines with Heterogeneous Servers," 

Operations Research ^ 8, 1960: 504-511. 

5. Koenigsberg, E. "Cyclic Queues," Operations Research . 9, 1958: 

22-35. 

6. Krishnamoorthi, B. "On Poisson Queue with Two Heterogeneous 

Servers," Operations Research , 11, 1963: 321-330. 

7. Makino, Toji. "On a Scheduling Problem," Journal of the Operations 

Research Society of Japan , 8. 1965: 32-44. 

8. McCullough, David U. "An Application of Queueing Theory to the 

Operation of Replenishment at Sea." United States Naval Post¬ 
graduate School, Monterey, October, 1966. 

9. Saaty, Thomas L. "Time Dependent Solution of the Many-Server 

Poisson Queue," Operations Research , 8, I960: 755-771. 

10. Smith, Richard D. and Dudek, Richard A. "A General Algorithm for 

Solution of the N-Job, M-Machine Sequencing Problem of the Flow 
Shop," Operations Research , 15, 1967: 71-81. 

11. Waggoner, Mark H. "An Analytical Model for Application to the 

Operation of Replenishment at Sea." United States Naval Post- 
Graduate School, Monterey, December, 1967. 


35 

























APPENDIX 


A. Conditional Laplace Transform for Steps 

X and Y are independent exponentially distributed random variables 
with parameters X and y respectively. Then 


P(X > Y) = If f(x) g(y) dy dx 
x>y 

00 X 

= / f (x) /. p e dy dx 

0 0 


= / Xe"^^ [1 - e-^^] dx 

0 


(Al) 


i 


P(X > Y) = 


P (X ^ 11X > Y) = 


P(Y < X 4 t) “ 


P 


X + ]i 


P(Y < X < t) 

P(X > Y) 


If f(x) g(y) dy dx 
y<x<t 

t * ■ _ 

/ f(x) / pe dy dx 

0 0 


/ Xe [1 - e"^*i“ dx 
0 


(A2) P(Y < X < t) = '■ -e'^^ + -r~Z — 

^ ^‘ 4 "X X*t*]i 

Combining the results of (Al) and (A2) yields 
(A3) P(X i t|x > Y) * 1 - e"^^ 


36 









Taking the derivative of this expression yields 




^x|x>Y*'^^ 


A(A + u) p -At -(A+y)t, 
-; Le - e J 


The conditional Laplace transform is then 


E[e ®^|X > Y] 


A (A + |i ) 

y 


1 

A + s 


1 

A + y + 


(A4) E[e ®^|X > Y] 


A(A + u) 

(A + s) (A + y + s) 
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